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ABSTRACT 
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The analysis of data from the Explorer 45 (S -A) electrostatic 
analyzer in the energy range 5 to 30 keV has provided some new results 
on the . *ng current ion composition. It has been well established that 
the storm time ring current has a decay time of several days during 
which the p.’rticle fluxes decrease nearly monotonically. In the past 
ring curreiil studies have assumed or stated that hydrogen was the dom- 
inant ion in the earth's ring current. By analyzing the measured ion 
fluxes during the several day storm recovery period and assuming that 
beside hydrogen other ions were present and that the decays were exponential 
in nature, we were able to establish three separate lifetimes for the 
ions. These fitted decay lifetimes are in excellent agreement with the 
expected charge exchange decay lifetimes for 0^, and He^ in the 
energy and L-value range of the data. This inference technique, thus, 
establishes the presence of measurable and appreciable quantities of 
oxygen and helium ions as well as protons in the storm-time ring current; 
we also find indications that He^^ may also be present under these same 
conditions. 



INTRODUCTION 


Since the early suggestion of Dossier and P arker (1959) that the 
charge exchange mechanism could be significant in the ring current 
decay process, this mechanism has varied between being in favor and out 
of favor in explaining ring current particle decay after magnetic storms. 

The long term decay times of several days for the measured fluxes were 
consistent with the available charge exchange lifetimes and, therefore, 
supported the charge exchange position ( Frank , 1967; Swisher and Fra 
1968). These favorable comparisons were based on the early computational 
work by Liemohn (1961), who specified both the mean lifetime, t^, for 
protons confined to the equatorial plane and the charge exchange lifetime, 
1 ^, for protons mirroring off the geomagnetic equator. More recently, 
other theoretical and observational aspects, such as the emergence of 
wave-particle theory (Cprnwal 1 et ^1 . , 1970) and the discrepancy pointed 
out by Lyons and Evans (1976) between the observed proton pitch angle 
distributions and those predicted by charge exchange decay weakened the 
chargp exchange position. It also became apparent that the overall measured 
decay lifetimes no longer agreed well with the more recently computed 
values ( Prol *^5 , 1973, Tinsley , 1976). 

Smi th et al . (1976) were among the first to use charge exchange 
lifetime values modified from those given by Liemohn (1961) for comparison 
to their satellite measurements. They employed both the neutral hydrogen 
density, n, of P rol ss (1973) and the proton charge exchange cross section, 

0 , of McUure (1966). L yons and Fvans (1976) used the equatorial charge 
exchange lifetimes given by Tinsley (1976), whose values were also a 
modification due u a better determination of the neutral hydrogen density 
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and charge exchange cross sections. Additional oroqress was made by 
S...i th and Bewtr a (1976) and independently by Cowle y (1977) in recomputing 
the relationship between the lifetimes of equatorially mirroring particles 
and those particles mirroring off the geomagnetic equator. Their results 
showed that the otf-90° pitch angle ion fluxes would not charge exchange 
decay as rapidly as was previously thought. In the review of the recent 
work associated with the determination of the charge exchange lifetimes 
for various ions in the magnetosphere. Smith and Bew tra (1977) provide 
the updated and corrected lifetime values especially through the ring 
current region. While this area of space research has been experiencing 
iome recent controversy over the various mechanisms, there is no question 
but that our understanding of the ring current source and loss mechanisms 
is progressing in a significant manner. 

In this paper we examine the consequences of making the assumption 
that, besides protons, ions such as helium and oxygen can repr.’sent a 
major fraction of the storm-time ring current particle flux. All of 
the improvements and new measurements in determining the charge exchange 
lifetimes (Stnith and Bewtra. 1977) are used. The charge exchange decay 
mechanism then is used to infer the ion comnosition in the ring current 
during the recovery phase of geomagnetic storms. 
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Background 


Until the recent work by Lyons and Evans (1976) and Tinsley (1976), 
in which they suggest He^ as a significant ion, ring current studies have 
assumed or stated that hydrogen was the dominant ion composing the ring 
current ( Berko et al . , 1975; and others before that). However, 0^ ions 
have been shown to be the dominant ion species under geomagnetic storm- 
time conditions at low altitudes ( Shelley et al . , 1972; Sharp pt al . , 

1976a, b). Shelley et al . (1976) report new observations of 0^ and fluxes 
streaming up the field lines from the ionosphere and Johnson et al . (1977) 
report observations of H^, He^ and 0^ in the ring current L-values at 
altitudes of 5000-7000 km. Indirect evidence of 0^ at higher altitudes 
has been offered by Mcllwain (1976) from observations at synchronous 
orbit on ATS-6 and by Frank et al . (1977) from IMP-7 at much larger al- 
titudes near the geomagnetic tail. 

In this paper we present the analysis of data from the Explorer 45 
(S -A) electrostatic analyzer. Ihis analyzer system measured the total 
ion flux in the energy per unit charge range of 1 to 30 keV. ^njjrh. ALA-I- 
(1976) used measurements from this same instrument in a previous analysis. 
Assuming that the instrument was responding primarily to protons, the 
observed flux decay rates during the February 24, 1972, magnetic stonn 
recovery generally agreed with tlie predicted charge exchange decay rates 
available at that time. There were, however, several difficulties. 

First, the observed decay rates, while in good agreement with those 
predicted by Li emoh n (1961), were longer than the newer predicted proton 
decay rales given by Prttl ss (1973). Second, the ordering of the data at 
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L = 3.5 was not as good as it was at the higher L-values of 4.25 and 5.0. 
Third, the off-90° pitch angle particles did not decay as rapidly as pro- 
ton charge exchange would have predicted 1961). 

Based on the work by Smith and Bewt ra (1976) regarding the predicted 
pitch angle dependence for the charge exchange decay mechanism and on 
the revised charge exchange lifetime values reviewed by Sui th and_Bewtra 
(1977), the Explorer 45 measurements have been re-analyzed. The 
following new procedures were employed: 

1. Only those particles having the same equatorial pitch angle were 

considered within a given data subset. Even with a near geographic 

equatorial orbit (3° inclination), the Exolorer 45 satellite does 

vary in magnetic latitude from orbit to on-it due to the 11° offset 

of the earth’s magnetic field axis (see Figure 7, Sini th et al . . 1976). 

This meant that for those of f-geomagnetic-equator orbits, fluxes 

measured at 90° local pitch angle did not mirror at X = 0°. In 

order to have a measurement every orbit (7.8 hrs.) at the same L-value 

and local time, X = 14° was the smallest mirror latitude which could 
m 

be used. A dipole geomagnetic field is assumed for all the calcula- 
tions in this paper. 

2. For the normal sampling rate of the electrostatic analyzer and the 
satellite spin period, two consecutive readouts corresponded roughly 

to pitch angles about 12° apart. These two readouts were linearly inter- 
polated to obtain the particle flux at the required pitch angle. Further, 
to reduce the statistical uncertainties, this process was repeated for 
the next such set of readouts for the same energy (usually a minute later) 
and the two calculated fluxes were averaged to obtain the flux at the 
specified L-value. 
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3. The neutral hydrogen density was determined individually for each of 
the geomagnetic storms considered. The geophysical conditions (solar 
activity, magnetic activity, season, etc.) for the specific time 
periods were used in this process. Smith and Bewtra (1977) provide 
the detailed description of this procedure. 

4. It was assumed that more than one ion species could be present in 
measurable quantities in the ring current region after these geomag- 
netic storms. The measured flux decay data were fitted with a double 
or triple exponential decay function of the form: 

f » E a.e'^^/^i^ (1) 

i ’ 

where the a^'s are tne fitted variables representing the relative 
quantities of each of the ion species, t is the time after the decay 
process begins, and the t^'s are the fitted variables representing 
the decay lifetime for each of the assumed ion species. In the 
functional fitting process the a^'s and t^'s were all free variables 
with no assumption made as to which specific ions were involved. 

5. The entire ring current recovery period was considered from the max- 
imum of the main phase to the time at which either the AE index or 
Dst was again enhanced by a disturbance. This will be discussed 
later in more detail. 
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FLUX DECAY OBSERVATIONS 


Magnetic Storm Conditions 

Ion measurements from the Explorer 45 (S -A) electrostatic analyzer 
made during three geomagnetic storms are presented in this section. Two 
of the storms, one beginning on February 24, 1972, and the other beginning 
on May 15, 1972, were moderate storms with maximum negative Dst (M. Sugiura 
and 0. J. Poros, private communication, 1973) of around 75y (Figures la 
and lb). Cah ill a nd Lee (1975) have described some of the Explorer 45 
magnetic field observations during the February storm. The third storm 
is the large, well known December 17, 1971, magnetic storm ( Cahi 1 1 , 1973) 
with a maximum negative Dst of l/ly (Figure Ic). The AE index provides 
a quantitative measure of the auroral zone magnetic activity produced by 
enhanced ionospheric currents flowing below and within the auroral oval 
(All en et a1 . , 1975). The trajectories of the satellite orbit for the 
three magnetic storms are shown in Figure 2. The magnetic local time 
(MLT) of apogee varied from 2100 MLT during the December storm to 1700 
MLT during the May storm. 

The measured ion fluxes at an energy of 9.2 keV and an L-value of 
4.25 outbound are shown in Figure 3 for the three storms for five or 
six day periods, which include the storm commencements, main phase maxima 
and the long recovery phases. Only those particles which mirrored at 
the same geomagnetic latitude of = 14° were selected. The equatorial 
* pitch angle corresponding to this mirror latitude is about 61.3°. 

In the February and May storms the particle fluxes at this energy 
increased nearly an order of magnitude after the sudden commencement and 
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did so within the 8 hours between one orbit and the next. Within this 
resolution the naxlmum observed flux was reached nearly simultaneously 
with the maximum of the main phase of the storm. The particle *‘lux then 
decayed monotonlcally over the next several days until new magnetic activity 
occurred, at which time the particle flux again Increased. In the December 
storm a similar pattern occurred with the exception of the Initial period 
after the main phase maximum. The extended period of "storm-related 
activity" (shown In Figure 3) was a result of large substorms occurring 
In conjunction with the principal main phase storm on December 17. This 
substorm activity can be identified from the AE index In Figure 1. It Is 
Interesting to note, however, that there was observed In the "storm-related 
activity" period a general reduction In particle flux for about a day and 
a half before a flux Increase occurred again on December 19. The AE Index 
(Figure 1) increased on this day to over 700 y. A period of "new activity" 
occurred on December 22 as Is evidenced by a flux increase (Figure 3) and 
by the AE and Dst activity (Figure 1) during thi: day. 

The following sequence of events is evident for these storms: 1) A 

period of principal storm-related activity occurred on a variable time 
scale depending on the specific geomagnetic storm, during which the ring 
current particle flux increased significantly and during which subsequent 
substorm injections may have kept the observed flux levels from decreasing 
in an unimpeded way; 2) next a decay period occurred during which new 
particle injections appeared to be negligible; 3) finally, new particle 
injections terminated the relatively quiet decay period. In the follow- 
ing, we shall be considering primarily the "decay period" as indicated in 
Figure 3. 
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Application of Charge Exchange Decay 
The improved charge exchange lifetimes given jy Smith an d Bewjra 
(1977) were used to determine the specific charge exchange lifetimes. 

T, for various ions at L-values and energies of interest. These values 
reflect the best deteniiination of the cross sections, neutral hydrogen 
density and mirror latitude dependence which are available as well as 
inputs of the specific geophysical conditions during the Explorer 45 
observations. These inputs are given in Table 1 (items 1-8) for the 
three magnetic storm periods. The detailed definitions of the density and 
temperature at the exobase, and T^, respectively, and the satellite 

critical altitude, are reviewed also by Smith and bewt ra (1977). Under 

the stated conditions in Table 1 (items 1-9) the neutral hydrogen densities 

n, were computed for the three storms and are given in item 10. 

For the ions H*, 0^, He^ and He^^ Figure 4 shows the normalized 

flux decay expected to be observed assuming only charge exchange decay 
and using the February , 1972, magnetic storm conditions. Ion energies of 
6.0, 9.2, 13.5 and 25.6 keV, and L-values of 3.5, 4.25 and 5.0 are con- 
sidered for the ions which are assumed to be mirroring at = 14°. In this 
paper we use the mass 4 isotope of the helium ions. For the singly charged 
ions, as can be seen from Figure 4, the<^e are significant and potentially ob- 
servable differences in the expected fluxes after 2-3 days of decay. In the 
electrostatic analyzer doubly charged ions, such as He , will be detected 
and recorded at energies which are half their actual energies. Therefore, 
to display the appropriate He^^ separation the cnarge exchange lifetimes 
used are for energies twice those shown. As can be seen in Figure 4 the He 
lifetimes are not different enough from H at all of these energies to always 
permit a separation of these two ions. 
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TABLE 1 : STORM PARAMETERS FOR CHARGE EXCHANGE D*'CAY 



Parameter 

Dec. 17, 1971 

Feb. 24, 1972 

May 15, 1972 

1. 

Ap 

11-19 (in 2 y) 

16-33 (in 2 y) 

17-21 (in 2 y) 

2. 

flO 7 flux 

133.8-143.8 

(10’22w/m^Hz) 

152.3-167.7 

(10'22y/m2„^) 

153.0-161.1 

llO-22w/m^Hz) 

3. 

81 day mean of f-|Q 7 

117.8 

128.4 

126.0 

4. 

Local time 

Assume symmetric 

Assume symmetric 

Assume symmetric 

5. 

Geographic Latitude 

Assume equatorial 

Assume equatorial 

Assume equatorial 

6. 

"c 

8.0 X 10^/cm^ 

6.5 X loVcm^ 

7.0 X 10^/cm^ 

7. 

Tc 

970°K 

1024°K 

1006°K 

8. 

"sc 

2.5 R^ 

2.5 R^ 

2.5 R^ 

9. 

L-value 

4.25 

4.25 

4.25 

TO. 

Density, n 

204/cm^ 

219/cm^ 

203/cm^ 



The charge exchange lifetimes of the 0*. and He Ions at these 

energies are different enough to theoretically allow Identification of 

different Ion species by measuring the decay rates of the Ions. The 

data shown In figure 3 during the "decay period" were fitted to a double 

exponential decay function (eg. 1). The a^'s (Ion fluxes) and t^'s 

(lifetimes) were free variables and determined by the oest fit to the data. 

♦ ♦ 

The charge exchange lifetimes of H and 0 for the same conditions as under 
which the data were taken agree remarkably well to fitted lifetimes and x^ 
In each of the three storms. For example. In the February storm x was 
3.27 X 10** sec and x was 1.63 x 10^ sec, whereas x(H ) Is 3.64 x 10 ser. 

4. 5 

and x(0 ) Is 1.78 x 10 sec. It appears, therefore, due to the close 
agreement of the charge exchange lifetime values with those fitted lifetime 
determined from '« actual measurements, that one can Infer the presence of 
significant quantities of 0^ In the storm-time ring current. This occurs 
not just In an Isolated case, but In each of the three storms we have anal- 
yzed. This inference technique was applied to the data at energies above and 
below the 9.2 keV channel described above as well as for L-values above and 
belr>< the L * 4.25 case also described above. 

Energy and L-Value Dependence 

The electrostatic analyzer on Explorer 45 has an energy range of 0.7 
to 30.2 keV. For convenience and to avoid consideration of coulomb scatter- 
ing effects which became more important below 1 to 2 keV we have considered 
four energy detector channels in the higher end of this energy range. The 
energy band width and center energy of these channels are given In Table 2. 
Figure 5 shows the measured fluxes at L = 4.25 in the four energy bands. 
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TABLE 2 : ENERGY BANDS 


Center Energy 
(keV) 


^min 

(keV) 

E 

max 

(keV) 

b.O 


5.25 

7.13 

9.2 


7.97 

10.81 

13.5 


11.74 

15.93 

25.6 


22.29 

30.23 
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6.0, 9.2, 13.5 and 25.6 keV for the February, 1972, geomagnetic storm. 

The data at 9.2 keV are identical to that shown in Figure 3. The same 
fitting procedure as for the 9.2 keV case described earlier was employed 
for the other three energies. The fitted values of t and t agree well 

I 2 

with the charge exchange lifetimes for and 0^. The solid curves 
through the data were determined by using the predicted charge exchange 
lifetime values for and 0^ with the fitted a and u values and normal - 

I 2 

ized to the first data point. Of course, curves derived using the fitted 

T and T values show even better agreement than these. 

1 2 

Since the Explorer 45 had an apogee slightly above L * 5 (Figure 2), 
we have selected three L-values (5.0, 4.25 and 3.5) which span the storm 
time ring current region as much as possible and are equally spaced. At 
L = 5, due to the altitude increase, the tjeutral hydrogen density, n, de- 

3 

creased and was computed to be 135/cm for the February storm. As described 

above, the data for each energy were fitted to a double exponential decay 

function (eq. 1) and the ion lifetimes, t^'s, and quantities, a^'s were 

determined by the best fit to the data. A comparison between the fitted 

i^'s and the charge exchange lifetimes of and 0^ for these same conditions 

showed remarkable agreement at L = 5. For example, at 13.5 keV = 5.70 x 

10^ sec and x = 2.26 x 10^ sec, whereas x(H^) is 5.88 x 10^ sec and t(O^) 

2 

is 2.31 x 10^ sec. 

The situation at L = 3.5 is somewhat different. The measured fluxes 
in the four energy bands at this L-value are shown in Figure 6. As in 
all cases before, the data for each energy were fitted to a double ex- 
ponential decay function. However, contrary to the previous cases, the 
chi-square of the fits were not nearly as good as before. For the 6.0 
keV and 9.2 keV cases it was found that better fits, based on chi-squared. 
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were obtained if the first data point (from February 24. 1972) was not 

included in the double exponential decay fit. The t and t values for 

4 S 

9.2 keV were found to be 9.34 x 10 sec and 7.90 x 10 sec which are in 

f > S 

close agreement to the 0 and He lifetimes of 1.09 x 10 sec and b.49 x 

B * 

10 sec, respectively. The situation is interpreted that the H decay 

at L * 3.5 during the first 8 hours wa*. substantial and thereafter only 

0^ and He^ existed in "observable" quantities. As can be seen in Figure 

4 at b.O keV and at L = 3.5 the drop in the normalized flux for during 

the first 8 hours under charge exchange decay is from 1.0 to 0.25, quite 

a substantial drop. 

The manner in which the above double exponential decay turves for 
b.O and 9.2 keV deviated from the measured fluxes lead us to consider 
using a triple exponential decay function of th*' form 


f - 


a P 1 

i 


It 




It e 


(t/i ) 

3 


(?) 


This implies that there could be thi'ee ions present in measurable 

quantities (if a / 0) with the i.'s sufficiently different from each 
^ 1 

other that the fitting procedure is sensitive to the difference. Amuzingly, 

5 

for 25. b keV the fitted value of i , 2.14 x 10 sec, is very close to the 

3 

predicted charge exchange lifetime* value He* ions would have, T(He*) = 

5 

2.40 X 10 sec, at this L-value and geophysical conditions where the 

3 

neutral hydrogen density, n, is 403/cm . The other two lifetimes, t and 

1 , were comparable to the charge exchange lifetime of H* and 0*, respec- 
2 

tively. A triple exponential decay was also used for the data at 13.5 keV 
and also yielded lifetimes close to the charge exchange decay lifetimes 
of H*, 0*, and He*. Apparently, with :he higher neutral hydrogen density 
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at this lower L-value the lifetimes of and are short enough such 
that they decay from their enhanced intensity to a level at which the 
He* intensity becomes "observable" through this technique. 

Kitrcr Latitude Dependence 

Thus far we have presented the Expli.rer 45 measurements of those 
particles which mirror near the geomagnetic equator. The relationship 
between the charge exchange lifetimes of particles mirroring off the 
geomagnetic equator to the lifetime of particles confined to the equator 
is ■ Tg cos"^ where j » 3.5 ± 0.2 v^^ith and Bewtf'a, 1976). For 
this February, 1972, magnetic storm at the L-values we are considering. 
Smith and Bewtra (1976) showed that the cosine dependence is more pre- 
cisely given by cos^'^X^. tor lifetimes, are reduced by 

only 10% from the equatorial values. However, for = 35“, for example, 
the lifetimes are shortened by nearly a factor of 2 from the equatorial 
values. The equatorial pitch angle for particles mirroring at 35° is 
about 28°. 

Figure 7 shows the measured ion fluxes at X^^ = 35° for the 9.2 keV 
energy channel at L » 4.25. A triple exponential decay function (eq. 2) 
was used to fit the data. The fitted values of t , t , and i are 2.81 

1 2 3 

4 5 5 

X 10 sec, 1.25 X 10 sec and 5.69 x 10 sec, respectively. These values 
are to be compared tc the H*, 0* arid He* charge exchange lifetimes, which 
are, 2.05 x 10*^ sec, 1.13 x 10^ sec and 6.69 x 10^ sec, respectively. 

The solid curve in Figure 7 is determined by using the predicted charge 
exchange lifetime values for the three ions and the fitted values for the 
ion quantities, a^'s. If the fitted t^'s were used the solid curve would 
again be in better agreement. 
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To show graphically the inferred contribution of He^ at these low 
equatorial pitch angles, the expected flux level of the combination of 
and 0*^ is shown as a dashed curve in Figure 7. This curve was de- 
termined by using the fitted a and a values, and the t(H*) and i{0*) 

J 2 

values. This show, that and 0 can not be the only two ions being 
measured. 
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OXYGEN AND HELIUM IONS IN RING CUnRENT 


The results of our analyses of the February 24, 1972, geomagnetic 
storm cleo»^ly infer that not only are oxygen ions abundant in the storm- 
time ring current, but helium ions arc also "observable". The helium 
ions in the energy range we are considering are most obvious at the 
lower altitude. We saw this for the ions at L ■ 3.5, the lowest L-value 
we considered as well as for the ions which were mirroring quite well 
down the field line (A^ ■ 35“). This is principally due to the fact 
that the charge exchange lifetimes are shorter for ions flowing through 
a larger density of neutral hydrogen such as occurs at the lower L-values 
and lower mirror latitudes, and the fitting technique for two or three 
ions is more sensitive to the shorter lifetimes since we have a limited 
number of measurements during a finite decay period (Figure 3). 

One is now confronted by the possibility that the double exponential 
decay fits, as good as they are at L « 4.25 and L = 5.0, do not yield 
the full informational content of the data. Therefore, we have fitted 
the triple exponential decay function to the data to provide consistency 
in our approach, realizing that the helium ions may be negligible in 
some of the cases. The three ion fit was applied to the February and 
May data (Figure 3) at L » 4.25 and fOi‘ an energy of 9.2 keV. We at- 
tempted to fit the December 17, 1971 storm at the same L-value and energy, 
but the "decay period" was too short to get a meaningful fit with the 
third ion. A minimum of three days is needed to get He^ fits under 
these conditions. The fitted x^'s and a^'s were used to determine the 
solid curves through the data in Figure 3 and are given in Table 3 for the 
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three storms. The lifetimes are also given for 0* and He^ at this 
energy and for the appropriate geophysical conditions given in Table 1. 

In all three stonns the fitted i/s agree well with the 0^ and He^ 
lifetimes. The value' generally are larger than the predicted values 
for and this difference will be addressed below in reference to a 
possible contribution. The other remarkable result, beside the 
agreement with the three ion species, is the general consistency of the 
fitted T.'s for the three stonns. 

The values represent the intensity of each of the ion species 
at t = 0, the time associated with the beginning of the fitting. This 
is approximately the beginning of the "decay period" as shown in Figure 
3. The t = 0 poiivc for the February and May storms is slightly delayed 
from the main phase maximum point. In the December storm, however, the 
t = 0 point is significantly shifted with respect to the main phase max- 
imum. The values of a , a and a for February and May are quite com- 

1 2 ) 

parable for the two storms and are approximately in the same pro- 
portions. In the December storm a is smaller than a . This inversion 

1 2 

of the relative quantities could be due to the fact that only a double 

exponential decay function was used in the fit, but it does not seem 

very likely. Instead what we may have observed is the charge exchange 

decay of ions (a ) before our t = 0 start to a level below 0^ ions 

1 

(a^). This would have occurred during the "storm-related activity" period 

in Figure 3. The values of a^, corresponding to He^ ions, are quite small 

at t = 0 compared to a and a . 

1 2 

The data at L = 5 are presented in Figure 8 for completeness of the 
data set for the February storm and the data are in the same format and 
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TABLE 3: ION LIFETIMES AND INTENSITIES (L ‘ 4.25; E - 9.2 keV)* 


Dec 17. 1971 Feb 24. 1972 May 15. 1972 


Fitted Lifetime t 

I 


Predicted Lifetime 


Fitted Lifetime i 

2 


Predicted Lifetime 0* 


Fitted Lifetime t 

3 


Predicted Lifetime He* 


2 _ 1 

a (cm stc sr keV) 

1 

2 _ > 

a (cm sec sr keV) 

2 

2 . 1 

a (cm sec sr keV) 

3 


*For 

m 


4.04 X 10^ sec 4.34 x 
3.91 X 10^ sec 3.64 x 
2.00 X 10^ sec 2.32 x 
2.16 X 10^ sec 2.01 x 

1.10 X 

1.28 X 10^ sec 1.19 x 
1.39 X 10^ 1.85 X 
2.06 X 10^ 1.28 X 

5.12 X 


10^ sec 4.22 x 10^ sec 

10^ sec 3.93 x 10^ sec 

10^ sec 2.17 X 10^ sec 

10^ sec 2.17 X 10^ sec 

10^ sec 1.10 X 10^ sec 

10® sec 1.29 X 10® sec 

10® 2.04 X 10® 

10® 1.06 X 10® 

10^ 5.06 X 10^ 
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with the same selection criteria for the data at L » 4.2b (Figure 5) 
and L ■ 3.5 (Figure 6). The solid curves, however, are the results of 
the triple exponential decay functional fits and are determined by using 
the predicted charge exchange lifetimes values for H^, 0* and He^ together 
with the fitted values for o , a and a as was done in the previous 

t 2 i 

cases. These curves show excellent agreement with the data. 

Up to this point we have presented the fitted lifetimes in comparison 
with the predicted charge exchange decay lifetimes on a case by case basis, 
i.e., for each given energy anJ L-value. In order that one gets a clearer 
understanding of the relative correspondence of the fitted t^'s to the 
predicted values, the fitted t/s determined by the triple exponential 
decay functions are shown as a function of energy in Figures 9a-c for the 
three L-values, 3.5, 4.25 and 5.0, respectively. The solid curves from 
1 to 100 keV are the predicted charge exchange decay lifetimes for the 
February 1972 geomagnetic storm for ions mirroring at = 14°. The 
fitted lifetimes at the four energies, 6.0, 9.2, 13.5 and 25.6 keV, cluster 
about the predicted lifetimes of H^, 0^ and He^ for each of the three 
L-values. The fitted lifetimes corresponding to 0^ and He^ ions are 
well ordered in energy, show excellent agreement with the predicted life- 
time values and are well separated from each of the other two ions. The 
one value for these two ions with the largest deviation from the pre- 
dicted value is in the 13.5 ke»/ energy band at L = 5.0. We have analyzed 
these data and can find no apparent reason outside possible substorm in- 
jection effects for the deviation, which is not all that large except 

t 

in comparison with the excellent agreement at all the other values. 

As was mentioned earlier the and He^^ charge exchange lifetimes 


19 



(with the energy of the He** shifted a factor of two) are close to one 
another and are overlapping at about 11 keV. The He** charge exchange 
lifetimes are shown also In Figure 9a-c, but with the energy scale shifted 
by a factor of two to reflect the energy bands In the electrostatic anal- 
yzer In which the He** Ions would be measured. The results In general 
for show overall consistency with H* especially If the He** curve 
Is Ignored. In examining this more closely, however, one Is drawn to the 
value at 25.6 keV which Is nearly In the middle of the H* and He** curves 
especially at the two lower L-values. At L * 5 It Is between the two 

curves but closer to the H* curve. In addition, the value of t at 6.0 

1 

keV Is also between the two curves. This strongly suggests a mixture of 
H* and He** Ions. Due to the relative closeness In the predicted life- 
times the fitting technique by Itself could not resolve He** from H*, 
as we were able to do with 0* and He*. The general trend of the fitted 
T^'s, considering the three L-values together suggest that at L * 5 the 
agreement Is much closer to H* and at L = 3.5 the trend has shifted toward 
a better agreement with He**. Recall also, that for the He** Ions the 
energy of the measured Ions (If any) Is twice the energy of the detector 
center energy. Thus any mixture of H* and He** which may be present are 
really mixtures at two different energies (e.g. 25.6 keV for H* and 51.2 
keV for He**). One can therefore consider that the He** may be more 
abundant at the higher energies than It is at the lower energies. 

Conclusions 

Through the application of the inference technique which we have 
described in this paper, it Is possible that the Ion composition of the 
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ring current In the energy range of an electrostatic analyzer can be 
Inferred for moderate to large geomagnetic storms. Using this method 
we Infer that oxygen Ions In the range 5 to 30 keV are quite abundant 
In the storm time ring current and that In the recovery phase, one or 
two days after the main phase maximum, 0* Ions can have Intensities 
greater than Ions. We also Infer that He^ can be found In "measurable" 
quantities at these energies and L-values but that the Intensities of 
He* Ions are much smaller than those for H* and 0*. However, as the 
recovery phase continues the H* Ions, which are removed rapidly by charge 
exchange decay, become small In Intensity, and the 0* Ions after several days 
are reduced to Intensities comparable to the He* Intensities, 

It thus appears that the He* Ions may form most of what we consider to 
be the quiet-time ring current, and It Is the hydrogen Ions which are the 
Subject of the Initial storm Injection and most assuredly are associated 
with substorir Injections. 

Additionally, the possibility that He** could also be present In 
the ring current is evident In our summary results. While the nature of 
this Identification Is not nearly as strong as our Inference of 0* and 
He*, It presents a distinct possibility. 

These results make one realize that Infrequent measurements through 
the ring current, with say an Ion mass spectrometer, will not sufficient- 
ly identify the ion composition since the mixture of Ions Is constantly 
changing due to their decay properties and due to the repeated Injection 
of fresh particles into this region. It would, however, give a definitive 
"snapshot" of the Ion composition which could be analyzed In terms of the 
scenerlo we have presented In this paper. 
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SUMMARY 


The analysis of data from the Explorer 45 (S^-A) electrostatic an- 
alyzer In the energy range 5 to 30 keV has provided some new results on 
the ring current Ion composition. By analyzing the measured decay rates 
of the Ions after geomagnetic storms and assuming two or more Ions were 
present and that the decays were exponential In nature we were able to 
establish three separate lifetimes for the Ions. By using the recently 
revised charge exchange decay lifetimes ( Smith and Bewtra , 1977) and the 
dependence of these lifetimes on the mirror latitudes of these particles 
we determined for three geomagnetic storms the expected charge exchange 
lifetimes for 0^ and He^ In the energy and L-value range of the data. 

In this paper we showed the close agreement of the fitted decay lifetimes 
from the measurements to the charge exchange decay lifetimes to these 
three Ions and Inferred that 0^ and He^ existed In the storm-time 
ring currents we considered. We also found Indication that He*^ may also 
be present. The agreement Is consistent In energy (5 to 30 keV), L-value 
(3.5 to 5.0), at least two widely spaced mirror latitudes * 14“ and 
35“ ) and for threa geomagnetic storms. We have thus developed an Inference 
technique to measure the ion composition In the ring current using an 
electrostatic analyzer. 


22 



REFERENCES 


Allen, J. H., C. C. Abston, and L. D. Morris, Auroral electrojet magnetic 
activity indices, AE(ll), for 1972, World Data Center A Report UAG-45, 
1975. 

Berko, F. W., L. J. Cahill, Jr., and T. A. Fritz, Protons as prime contri- 
butors to the storm-time ring current, J. Geophys. Res ., 80, 3549- 
3552, 1975. 

Cahill, L. J., Jr., Magnetic storm inflation in the evening sector, 

J . Geophys. Res ., 78 , 4724-4730, 1973. 

Cahill, L. J., Jr., and Y. C. Lee, Development of four small magnetic 
storms in February, 1972, Planet. Space Sci . , 23, 1279-1292, 1975. 

Cornwall, J. M. , F. V. Coroniti, and R. M. Thorne, Turbulent loss of 
ring current protons, J . Geophys. Res ., 75, 4699-4709, 1970. 

Cowley, C. W. H., Pitch angle dependence of the charge-exchange lifetime 
of ring current ions. Planet. Space Sci . , 25 , 385-393, 1977. 

Dressier, A. J., and E. N. Parker, Hydromagnetic theory of geometric 
storms, J. Geophys. Res ., 64 , 2239, 1959. 

Frank, L. A., On the extraterrestrial ring current during geomagnetic 
storms, J. Geophys. Re s., 72 , 3753-3767, 1967. 

Frank, L. A., K. L. Ackerson, and D. M. Yeager, Observations of atomic 
oxygen (O"*^) in the earth's magnetotail, J. Geophys. Res ., 82, 129- 
134, 1977. 

Johnson, R. G., R. D. Sharp, and E. G. Shelley, Observations of ions of 
ionospheric origin in the storm-time ring current, Geophys . Res . 
Letters , 4, 403-406, 1977. 

Liemohn, H., The lifetimes of radiation belt protons with energies be- 
tween 1 keV and 1 MeV, J. Geophys. Res ., 66, 3593-3595, 1961. 

Lyons, L. R., and D. S. Evans, The inconsistency between proton charge 
exchange and the observed ring current decay, J. Geophys. Res ., 

§1, 6197-6200, 1976. 

McClure, G. W., Electron transfer in proton-hydrogen-atom collisions: 

2-117 keV, Phys. Rev ., 148 , 47-54, 1966. 

Mcllwain, C. E., Ionospheric ions at 6.6 earth radii, in E^, 1212, 

1977. 

Prolss, G. W., Decay of the magnetic storm ring current by the charge- 
exchange mechanism. Pl anet. Space Sci . , 21 , 983-992, 1973. 


23 



Sharp, R. 0., R. G. Johnson, and E. G. Shelley, The morpho’rgy of ener- 
getic 0* ions during two magnetic storms: Temporal variations, J . 

Geophys. Res., 81 , 3283-3291, 1976a. 

Sharp, R. D., R. G. Johnson, and E. G. Shelley, The morphology of ener- 
getic 0* Ions during two magnetic storms: Latitudinal variations, 

J. Geophys. Res ., 6 ]_, 3292-3298, 1976b. 

Shelley, E. G. , R. G. Johnson, and R. D. Sharp, Satellite observations 
of energetic heavy Ions during a geomagnetic storm, J. Ge ophys. 

Res. , 21 , 6104-6110, 1972. 

Shelley, E. G., R. 0. Sharp, and R. G. Johnson, Satellite observations 
of an ionospheric acceleration mechanism, Geophys. k es. L e tters . 

3, 654-656, 1976. 

Smith, P. H., and N. K. Bewtra, Dependence of the charge exchange life- 
times on mirror latitude, Ge ophys. Res. Letters . 3, 689-692, 1976. 

Smith, P. H., and N. K. Bewtra, Charge exchange lifetimes for Ions In 
the magnetosphere, (GSFC Preprint X-620-77-209) , 1977, submitted 
for publication. 

Smith, P. H., R. A. Hoffman, and T. A. Fritz, Ring current proton decay 
by charge exchange, J. Geophys. Res ., 81 , 2701-2708, 1976. 

Swisher, R. L., and L. A. Frank, Lifetimes for low-energy protons In the 
outer radiation zone, J. Geophys. Res ., 73 . 5665-5672, 1968. 

Tinsley, B. A., Evidence that the recove-y phase ring current consists 
of helium ions, J. Geophys. Res ., 81 . 6193-6196, 1976. 


24 



FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Ost (y) and the AE Index (y) For six day periods covering 
three geomagnetic storms: December 17, 1971, February 24, 

1972, and May 15, 1972. 

Trajectories of the Explorer 45 (S^-A) satellite orbit in 
L versus MLT coordinates during the three geomagnetic storms. 

The trajectories are for the near-equatorial orbits during 
these periods. The data presented In this paper were taken 
during the outbound portion of the orbits. 

Ion fluxes measured by the Explorer 45 electrostatic analyzer 
at L » 4.25 with an energy of 9.2 keV for the three geomagnetic 
storms. The fluxes were of Ions mirroring at the geomagnetic 
latitude of 14®. This corresponds to an equatorial pitch angle 
of 61.3°. The sudden commencements (SC) and the times of the 
main phase maximum (MPM) are indicated. The solid curves 
through the data points are explained later In the text. 

Nornialized flux decay expected for the Ions H^, 0*, He^ and 
He^'*’ for four energies (6.0, 9.2, 13.5 and 25.6 keV) at three 
L-values (3.5, 4.25, 5.0) assuming only charge exchange decay 
mechanism of Ions mirroring at * 14° and using the February 
24, 1972 magnetic storm c^ijidltlons to determine the neutral 
hydrogen density. The He^^ decay curves are for Ions at twice 
the energy Indicated. 

Ion fluxes measured by the Explorer 45 electrostatic analyzer 
at L = 4.25 in four energy channels (6.0, 9.2, 13.5 and 25.6 
keV). The fluxes were of ions mirroring at the geomagnetic 
latitude of 14°. The solid curves, described In detail in the 
text, were determined by assuming that the data represents a 
mixture of particle flux of two ions. 

Same as Figure 5 except now for L = 3.5. The solid curves, 
described in detail in the text, were determined for 6.0 and 
9.2 keV by assumluy that the data represents a mixture of 
particle flux of t’/o ions. The curves through the 13.5 and 
25.6 keV data are determined by assuming a mixture of three 
ions. 

Ion fluxes at L = 4.25 but for those 9.2 keV-particles mirror- 
ing at Xn, = 35°. The solid curve was determined by assuming 
a mixture of three ions and is described in the text. The 
dashed curve Is the normalized charge exchange decay curve 
under the stated conditions for a combination of and O’*" ions. 
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